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Abstract: I discuss some aspects of the recent test of frame-dragging performed by me by exploiting the Root- Mean- 

Square (RMS) orbit overlap differences of the out-of-plane component N of the orbit of the Mars Global 
Surveyor (MGS) spacecraft in the gravitational field of Mars. A linear fit of the full time series of the entire 
MGS data (4 February 1999-14 January 2005) yields a normalized slope 1.03 ± 0.41 (with 95% confidence 
bounds). Other linear fits to different data sets confirm the agreement with general relativity. The huge 
systematic effects induced by the mismodeling in the martian gravitational field claimed by some authors 
are absent in the MGS out-of-plane record. The non-gravitational forces affect at the same level of the 
gravitomagnetic one the in-plane orbital components of MGS, not the out-of-plane one. Moreover, they 
experience high-frequency variations which does not matter in the present case in which secular effects are 
relevant. 
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1. Introduction 

I proposed in Refs. [1, 2] an interpretation of the time series of the RMS orbit overlap differences [3] of the 
out-of-plane part N of the orbit of the Martian polar artificial satellite Mars Global Surveyor (MGS) over a time 
span AP of about 5 years (14 November 1999-14 January 2005 in Ref. [2]) in terms of the general relativistic 
gravitomagnetic Lense-Thirring effect 1 [6] . It turned out that the average of such a time series over AP, normalized 
to the predicted Lense-Thirring out-of-plane mean shift over the same time span, is fi — 1.0018 ± 0.0053. 
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It consists of a small precession of the orbit of a test particle freely moving around a central rotating mass; 
the general relativistic component of the gravitational field due to the mass currents generated by the rotation of 
body is named "gravitomagnetic" by analogy with the magnetic field induced by electric currents in the Maxwellian 
electromagnetism [4J. For some recent, controversial tests performed in the gravitational field of the Earth with 
the LAGEOS satellites, see Ref. [5] and references therein. 
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Our interpretation has recently been questioned by Krogh in Ref. [7]. The remarks concerning the analysis 
presented in Refs. [1, 2] mainly deal with I) The observable used: in Refs. [1, 2] I would have misinterpreted the 
MGS data II) The confrontation between the prediction of the gravitomagnetic Lense-Thirring shift and the data 
over the chosen time span AP: in Refs. [1, 2] I would have incorrectly compared the 1.6 m value of the out-of- 
plane average orbit error released by Konopliv et al. in Ref. [3] for the entire MGS data set to the Lense-Thirring 
shift calculated for a shorter time interval AP III) The data set used: in Refs. [1, 2] I discarded some of the initial 
months of the MGS data set IV) The bias-neglected by me in Refs. [1, 2]-due to the multipolar expansion of the 
Newtonian part of the martian gravity field, as pointed out in the unpublished Refs. [8, 9] quoted by Krogh in 
Ref. [7] V) The impact of the atmospheric drag, neglected by me in Refs. [1, 2] 

Below I present my reply which, basically, consists of the following points. As further, independent tests, here I 
present various linear fits to different data sets including, among others, the full time series of the entire MGS data 
(4 February 1999-14 January 2005) as well; the predictions of general relativity turn out to be always confirmed. 
The analytical calculation of the competing aliasing effects due to both the gravitational and non-gravitational 
perturbations, which affect the in-plane orbital components of MGS, do not show up in the real data. Moreover, 
the non-conservative forces, whose steadily refined modeling mainly improved the in-plane orbital components of 
MGS, not the normal one, exhibit high-frequency, non-cumulative in time variations. 



Before going into the details of my reply, I will extend to the MGS-Mars system a standard argument used to 
evaluate the possibility that general relativity has to be taken into account in describing the motion of test bodies 
in a given gravitational field [10]. It consists of considering the relativistic characteristic length of the problem 
at hand, like the Schwarzschild radius of the body acting as source of the gravitational field, and comparing it 
to the accuracy with which the orbit of a test particle orbiting it can be determined. In the present case, the 
characteristic gravitomagnetic length of a rotating body of mass M and angular momentum S is [11] 

<• = m> (1) 

where c is the speed of light in vacuum; l g is also one of the two parameters with dimension of lengths entering 
the Kerr metric [12] which describes the exterior field of a rotating black hole. Since the angular momentum of 
Mars can be evaluated as 



from the latest spacecraft-based determinations of the areophysical parameters quoted in Ref. [3], it turns out 



2. My arguments 



2.1. A consideration of general nature 



S M = (1.92 ± 0.01) x 10 32 kg m 2 s 



-l 



(2) 




cM M 



= 1.0 m. 



(3) 
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Such a value has to be compared with the present-day accuracy in determining the orbit of MGS; it can be 
evaluated it is about 0.15 m [3] in the radial direction, not affected by the gravitomagnetic force itself, as I 
will show later. Thus, it makes sense to investigate the possibility of measuring the Lense-Thirring effect in the 
gravitational field of Mars with MGS in such a way that a possible positive outcome should not be regarded as 
unlikely. 

2.2. Reply to Krogh's points 

1) The entire MGS data set was subdivided by Konopliv et al. in Ref. [3] in 388 (not 442, as claimed 
by Krogh in Ref. [7]) smaller time intervals of data called arcs. For MGS, the lengths of the arcs vary 
from 4 to 6 days, so to cover many orbital revolutions (~ 2 h). For each arc, the spacecraft position and 
the velocity, among other things, were estimated and used as starting point for a numerical propagation 
of the satellite's motion by means of the dynamical models which, in the case of MGS, did not include 
the general relativistic gravitomagnetic force. Contiguous arcs were overlapped by an amount of just 2 
h, i.e. one orbital revolution, and the 2 RMS spacecraft position difference among the predicted positions 
propagated from the estimated ones in the previous arc and the estimated positions of the subsequent arc 
was computed. Since the arc overlaps cover just about one orbit, such RMS differences, in fact, account for 
any among measurement errors, random errors, systematic bias due to mismodeling/unmodeling dynamical 
forces yielding secular, i.e. averaged over one orbital revolution, effects, whatever their physical origin may 
be. Indeed, RMS of orbit solution overlaps are commonly used in satellite geodesy as useful and significant 
indicators of the overall orbit accuracy [13, 14]. Conversely, they are also used to gain information about 
systematic errors coming from inaccurate modeling of the forces acting on the spacecraft. For details see 
Refs. [13, 14]. Of course, such a technique is insensitive to short-period effects, i.e. having frequencies 
higher than the orbital one: only dynamical features of motion with time-scales equal to, or larger than one 
orbital period can be sensed by such orbit overlap differences. Moreover, the average orbit error (AiVdifr) 
of about 1.6 m does not refer to this or that particular arc overlap; instead, it comes from the mean of 
the entire set of RMS orbit overlap differences for the chosen time span AP and is well representative of 
those un-modelled/mis- modelled forces yielding effects which do not average out over AP, as it is just the 
case of the Lense-Thirring signal. Time- varying patterns exhibiting well-defined periodicities-including also 
measurement errors like, e.g., those related to the Earth-Mars geometry-are, instead, mainly averaged out 
yielding little or no contribution to the average orbit error. Incidentally, from the above discussion about 
the meaning of the average orbit error, it should be apparent that it does not make sense to look for the 

2 / acknowledge the use of wrong terminology in [1], as noted by Kroght in Ref. [7], although I feel that the 
complaint of plagiarism that seems to be raised by Krogh in Ref. [7] in his footnote 1, pag. 5710 against me [2] 
about such a matter sounds a bit excessive. 
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error of the error, as, instead, seemingly required by Krogh in Ref. [7] when he blames me [2] for not 
having included the uncertainty in (AN TCB ). Another criticism by Krogh in Ref. [7] is that the RMS overlap 
differences would be unable to specify any orbital precession. 

To reply to all such criticisms I decide to perform another, independent test of my hypothesis. First, by 
linearly fitting 3 the full time series of Ref. [3], after having rescaled the data points in order to shift the 
zero point of the time-series to the middle of the data span, I get a slope of —0.64 ± 0.26 m yr" 1 , (with 
95% confidence bounds), while the predicted Lense-Thirring MGS out-of-plane rate (customarily defined 
positive along the spacecraft's orbital angular momentum) amounts to 0.62 m yr -1 . The obtained minus 
sign is due to the fact that Konopliv et al. in Ref. [3] defined the normal direction to be positive in the 
opposite direction of the MGS orbital angular momentum (Konopliv 2007, private communication) . Should 
such a linear fit be used as indicator of the existence of the Lense-Thirring effect, its relativistic prediction 
would be fully confirmed within the experimental error; instead, the hypothesis of a null effect would be 
rejected at 2.4 sigma level. Then, I also repeat my procedure by fitting with a straight line the entire data 
set without full January 2001, mainly affected by likely measurement errors which, according to Krogh in 
Ref. [7], would mimic the Lense-Thirring effect, getting —0.61 ± 0.26 m yr . The removal of the entire 
year 2001, mainly affected by angular momentum wheel desaturation operations, yields —0.57 ± 0.28 m 
yr" 1 . Another linear fit to the time series after removing the last month (December 2004- January 2005) 
yields -0.62 ± 0.27 m yr" 1 . 

Such results reply to the criticisms II) and III) as well concerning AP, to which a large part of Ref. [7] is 
devoted. 

IV) Krogh in Ref. [7] quotes Ref. [8] in which analytical calculation about the corrupting impact of various 
physical parameters of Mars through the classical node precessions induced by the even zonal harmonic 
coefficients Ji of the multipolar expansion of the Newtonian part of the martian gravitational potential are 
presented. In particular, Sindoni et al. in Ref. [8] use the first five even zonals J2...J10 along with their 
associated errors from former global solutions for the Mars'gravity field, the uncertainty in the Mars' GM 
and in the MGS semimajor axis and inclination, plug them into analytical formulas for the classical secular 
node precessions and conclude that, since the resulting effect is tens of thousand times larger than the Lense- 
Thirring effect on MGS, this would be fatal for any attempt to detect the gravitomagnetic frame-dragging 
with such a spacecraft. 

The point is that such figures, as others which can be obtained from more accurate calculation, must 
ultimately be compared with the reality of the data, i.e. the RMS orbit overlap differences of MGS. 

3 Note that, since the plots in Fig. 3 of Ref. [3] are semi-logarithmic, one should not visually look for a straight 
line in them. 
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I, in fact, repeated such calculation by considering also the other even zonals up to J20 along with the 
latest errors of the MGS95J global solution and including the uncertainties in the Mars' radius as well. By 
summing, in a root-sum-square fashion, such terms I get a mean bias of 78.9 m d _1 in the out-of-plane 
MGS orbital component: by linearly summing them I get an upper bound of 111.6 m d _1 . Such figures 
clearly show how they are by far not representative of the real MGS orbit. Indeed, over a time span of 5 
years I would have an enormous mean shift as large as 144 ^ m (root-sum-square calculation) or 203 km 
(linear sum). Interestingly, even if the set of the RMS overlap differences of MGS were to be considered 
as representative of a single orbital arc 6 d long only, the conclusion would be the same: indeed, in this 
case, the total cross-track mean shift due to the martian gravitational potential would amount to 473.1 m 
(root-sum-square) or 669.6 m (linear sum). 

In regard to Ref. [9] , quoted in Ref. [7] as well, let us recall again that the RMS orbit overlap differences are 
just used to account, in general, for all the measuremnt/systematic errors giving an indication of the overall 
orbit accuracy [13, 14]. The important point is that they cancel out, by construction, errors, systematic or 
not, common to consecutive arcs— it would just be the case of a bias like that described by Felici in Ref. 
[9] — , while effects like the Lense-Thirring one, accumulating in time, are, instead, singled out [14]. 

V) In regard to the impact of the non-gravitational perturbations, Sindoni et al. in Ref. [8] yield a total 
un-modelled non-gravitational acceleration of ~ 10 -11 m s -2 which is the same order of magnitude of the 
Lense-Thirring acceleration induced by Mars on MGS. They neither present any detailed calculation of the 
effect of such an acceleration on the normal portion of the MGS orbit nor specify if such a magnitude refers 
to the out-of-plane component. However, some simple considerations can be easily traced: a hypothetic, 
generic perturbing out-of-plane force 6.7 times larger than the Lense-Thirring one and having the same 
time signature, i.e. linear in time, should induce a 10.8 m cross-track shift, on average, over the considered 
time span AP. Again, such a bias is neatly absent from the data. By the way, as clearly stated in Ref. [3], 
it is the along-track portion of the MGS orbit— left unaffected by the Lense-Thirring force— to be mainly 
perturbed by the non- gravitational forces: indeed, the along-track empirical accelerations fitted by Konopliv 
et al. in Ref. [3] amount just to w 10 -11 m s~ 2 , which shows that the guess by Sindoni et al. in Ref. [8] is 
somewhat correct, but it refers to the along-track component. 

Time-dependent, periodic signatures would, instead, be averaged out, provided that their characteristic 
time scales are relatively short, as it is just the case. Indeed, the non-conservative accelerations, which are 
especially active in the MGS in-plane orbital components as clearly stated in Refs. [3, 15], exhibit time- 
varying patterns over 12 hr [15] which, hypothetically mapped to the out-of-plane direction, are averaged 
out over multi-year time spans (and, incidentally, over 6 d as well). To be more definite, in regard to the 
issue of the impact of the atmospheric drag on the cross-track portion of the orbit of MGS, raised by Krogh 
in Ref. [7], let us note that it requires not only to consider the node Q, as apparently claimed by Krogh in 
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Ref. [7], but also the inclination i according to Ref. [16] 

AN = aJ(l+^) [ ( ^ + (siniAfi) 2 . (.1) 



According to, e.g., Milani et al. in Ref. [17], the perturbing acceleration Adrag due to the atmospheric drag 
can be cast into the form 

1 S 

^drag = --ZCv—pvv, (5) 

where S/M is the spacecraft cross sectional area (perpendicular to the velocity) divided by its mass, Co is 
the drag coefficient, p is the atmospheric density (assumed to be constant over one orbital revolution), v 
is the satellite velocity in a planetocentric, non-rotating frame of reference and Z is a corrective coefficient 
accounting for the fact that the atmosphere is not at rest, but rotates with angular velocity uja more or 
less rigidly with the planet; Z ~ 1 for polar orbits [17]. While the secular, i.e. averaged over one orbital 
period T, drag shift on the node vanishes, it is not so for the inclination: indeed, it turns out [17] 

{^ T ~J^p)^ + 0{el (6) 
\ n z a J n 

where n = \J GM/a 3 is the Keplerian mean motion. As a result, the orbital plane tends to approach the 
planet's equator; the terms in brackets is the ratio of the drag force to the Newtonian monopole. As usual 
in perturbation theory, o is meant as evaluated on the unperturbed reference ellipse. Thus, the out-of-plane 
drag shift is from (4) 

(AiVdrag ) « a <^p>. (7) 
In the following I will assume that uja ~ WMars ~ 

7.10 x 10" 5 s" 1 . Let us see what happens in the (unlikely) 



<AiV drag ) T ~ 1 x 10 -5 m. (8) 



But Adrag is not constant over time spans days or years long [18], so that such an effect is not a concern 
here. By the way, even if it was not so, by assuming a m 10% mismodeling in drag-which is, in fact, 
modeled by Konopliv et al. in Ref. [3]-(8), mapped onto about 5 yr, would give a w 0.7% uncertainty. 

Finally, Krogh in Ref. [7] remarks that decreasing in the averages of the RMS orbit overlaps occurred in 
view of constantly improved modeling [15, 19], but he does not recognize that the improved modeling of 
the non-gravitational forces acting on MGS introduced in Ref. [3] with respect to previous works [15, 19] in 
which the Lense-Thirring effect was not modelled as well, only affected in a relevant way just the along-track 
RMS overlap differences (a factor 10 better than in Ref. [15, 19]), not the normal ones (just a factor 2 
better than in Ref. [15, 19]). Moreover, if the relativistic signature was removed or not present at all so 
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that the determined out-of-plane RMS overlap differences were only (or mainly) due to other causes like 
mismodeling or unmodeling in the non-gravitational forces, it is difficult to understand why the along-track 
RMS overlap differences (middle panel of Figure 3 of Ref. [3]) have almost the same magnitude, since the 
along-track component of the MGS orbit is much more affected by the non-gravitational accelerations (e.g. 
the atmospheric drag) than the out-of-plane one. 
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